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Abstract. We present explicit expressions for the lens equation for a cored isothermal ellipsoidal gravitational 
lens as a a single real sixth-order algebraic equation in two approaches; 2-dimensional Cartesian coordinates and 
3-dimensional polar ones. We find a condition for physical solutions which correspond to at most five images. For 
a singular isothermal ellipsoid, the sixth-order equation is reduced to fourth-order one for which analytic solutions 
are well-known. Furthermore, we derive analytic criteria for determining the number of images for the singular 
lens, which give us simple expressions for the caustics and critical curves. The present formulation offers a useful 
way for studying galaxy lenses frequently modeled as isothermal ellipsoids. 
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1. Introduction 

Gravitational leasing due to a galaxy is important for 
probing mass distributions and determining cosmologi¬ 
cal parameters. Galaxy lenses are often modeled as cored 
isothermal ellipsoids. Although the ellipsoidal model is 
quite simple, it enables us to understand a number of 
physical properties of the galactic lens. Furthermore, it 
fits well with mass profiles implied by observations (For 
instance, Binney and Tremaine 1987). 

Until now, the lens equation for the cored isothermal 
ellipsoid has been solved numerically as a nonlinearly cou¬ 
pled system. For a binary gravitational lens, it has recently 
been shown that the lens equation is reduced to a single 
real fifth-order algebraic equation (Asada 2002). Also for 
a singular isothermal ellipsoidal lens, furthermore, the ap¬ 
parently coupled lens equations can be reduced to a single 
equation (Schneider et al. 1992), though explicit expres¬ 
sions were not given there. Along this course, we reex¬ 
amine the coupled lens equations for the cored isothermal 
ellipsoid. The main purpose of the present paper is to show 
that they are reduced to a single equation with a condition 
for physical solutions and to give analytic criteria for de¬ 
termining the number of images for isothermal ellipsoidal 
lenses. 
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2. Lens Equation for a Singular Isothermal 
Ellipsoid 

First, let us consider the singular isothermal ellipsoidal 
lens with ellipticity 0 < e < 1/5. A condition that the 
surface mass density projected onto the lens plane must 
be non-negative everywhere puts a constraint on the ellip¬ 
ticity as e < 1. A tighter constraint e < 1/5 comes from 
that the density contours must be convex, which is rea¬ 
sonable for an isolated relaxed system. The lens equation 
is expressed as 
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where (3 = (/3i,/32) and 6 = {01,02) denote the positions 
of the source and images, respectively. 

For simplicity, we introduce variables as cc = y/1 — e0i , 
y = v^l -I- €02, a = VI — e/3i and b = \/\ A- e/32, so that 
the lens equation can be rewritten as 
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2.1. Sources on the symmetry axes 


2.2. Off-axis sources 


There are two symmetry axes in the ellipsoid, a = 0 and 
6=0. We consider a source on the axis a = 0. In this 
case, we can find analytic solutions for the lens equation 
as follows. We can make a replacement as a 6 , x <-!• y 
and l + ee-i-l — eto obtain solutions for 6 = 0 . 

For 0 = 0, Eq. (|3|) becomes 


X 



^/aFT^ j 


= 0 , 


which means 
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a: = 0, 
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or 


In the case of a; = 0, Eq. (^) becomes 
6 -y = -(I + e)|^, ( 8 ) 

which is solved as 

r 6 -(l + e) for 6 <-(I + e), 
y = < 6 ± (1 + e) for —(1 + e) < 6 < 1 + e , (9) 

[ 6 + (1 + e) for 1 + e < 6 . 

Here, we should pay attention to the cases of 6 = ±(I + e), 
since they produce a solution at the singularity (a;, y) = 
( 0 , 0 ) of the potential. 

Next, we consider the case of Eq. (^. Let us investigate 
the two cases, e = 0 and e ^ 0 separately. For e = 0, Eq. 
means 6 = 0 , so that images become a ring as 

a ;2 + y 2 ^ 1 . ( 10 ) 


Below, we assume e 7 ^ 0. Eliminating \/x'^ y'^ from Eqs. 

and (|^), we obtain 
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Substituting this into y in Eq. (Q), we obtain 
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which has real solutions if and only if | 6 | < 2 e, 

X = ± ^ \/ 4e2 — 62. (13) 

Consequently, Eqs. (^ and ( |l^ show that four, two or one 
images occur for | 6 | < 2 e, 2 e < | 6 | < 1 + e or 1 + e < | 6 |, 
respectively. 

In the similar manner, we obtain the solutions for 6 = 
0. A point is that 2e < |a| < 1—e can hold only for e < 1/3. 


Here, we consider off-axis sources (a 7 ^ 0 and 6 7 ^ 0). 
In this case, Eqs. (^ and (Q) show a; 7 ^ 0 and y ^ 0. 
Eliminating \/x'^ y'^ from Eqs. (^) and (jj), we obtain 


(1 — e)bx 

^ (1 -f e)a — 2ex ’ 


(14) 


which determines y uniquely for any given x. For finite 
6 , Eq. (^) means that y is also finite, so that a; 7 ^ (1 -I- 
e)a/2e from Eq. &■ Substituting Eq. (0) into Eq. (^, 
we obtain the fourth-order polynomial for x as 


D{x) = [(a — a;)2 — (1 — e)2][(l-|-e)a — 2ea;]^ 

-|-(1 - e)26^(a - x)"^ 

= 0, (15) 


where we used a: 7 ^ 0 . 

The number of real roots for a fourth-order equation 
is discussed by the discriminant D 4 (e.g. van der Waerden 
1966), which becomes for Eq. ( 0 ) 

D 4 = -640^6^6^(1-e)^2 

X [(a2 -h 62 - 46^)3 -F 108a^b^e% (16) 


Namely, if 
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the number of real roots is either four or zero, which is de¬ 
termined as four by explicit solutions for on-axis sources. 
Otherwise, it is two. However, the number does not nec¬ 
essarily indicate that of images as shown below. 

Since a; 7 ^ 0 for off-axis sources, Eq. (i is rewritten as 

^ = - (18) 
^ \Jx^ y y2 

whose right-hand side is necessarily negative since 1 — 
6 > 0. As a result, we find that any solution of the lens 
equation must satisfy 


a — X 
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< 0 . 
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This implies that a; < 0 or a < a; for positive a, while 
a; < a or 0 < a: for negative a. It should be noted that Eq. 
( |l^ always holds in the limit of a ^ 0 . 

Let us investigate the number of roots in (a—x)/x > 0, 
namely an interval between 0 and a. We hnd out 


D{a) = - 0^(1 - ef < 0, 
D{Q)=a\l-e^f 
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If a and 6 satisfy 
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11(0) is not negative, so that D{x) = 0 has at least one 
root between 0 and a. Eor e < 1/3, Eq. (p^ implies 
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Fig. 1. Caustics for a singular isothermal ellipsoidal lens 
in the physical coordinates (/ 3 i,/ 32 )- For e = 0.15, the solid 
curve denotes the inner caustic, and the dashed one for 
the outer caustic. Sources locate at (0,0.1), (0,0.5) and 
(0,1.3), denoted by the circle, filled disk and square, re¬ 
spectively. 
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a^/4e^ -I- 6^/4e^ > 1, so that the left-hand side of Eq. © 
becomes positive. Hence, D 4 is negative, so that D{x) = 0 
has two roots. As a result, it has only one root in the in¬ 
terval, which means only one image appears. Unless Eq. 
( p^ ) holds, D{0) is negative, so that the number of roots 
for D{x) = 0 for (a — x)lx < 0, namely that of images are 
four or two, respectively for D 4 > 0 or < 0 . 

The inner and outer caustics (Eig. 1) are given by 
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The inner caustic given by Eq. ( ^^ 
is parametrized as 

a = 2 e cos^ t, 
b = 2 e sin^ t. 


(24) 

is an asteroid which 
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where t G [0, 27r). 

The critical curves on the lens plane correspond to the 
caustics on the source plane (Schneider et al. 1992). We in¬ 
troduce the polar coordinates as (cc, y) = (pcos^,psin^). 
By substituting Eqs. ( |^) and ( g 6 |) with t = into the 
lens equations (^ and^), we obtain the parametric rep¬ 
resentation of the critical curve as 
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Fig. 2. Critical curves for a singular isothermal ellipsoidal 
lens in the physical coordinates (0i, 02 )- For e = 0.15, the 
solid curve denotes the outer critical curve which corre¬ 
sponds to the inner caustic. The origin corresponds to the 
outer caustic. The images correspond to the sources in 
Fig. 1. 
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3. Lens Equation for a Cored Isothermal Ellipsoid 

Let us consider a cored isothermal ellipsoidal lens with the 
angular core radius c. The lens equation is expressed as 


a = a; 1 — 


b = yll- 


1 -e 


-I- y2 _|_ £.2 J 
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which are apparently similar to the set of Eqs. (H) and (^ . 
However, there are differences in their algebraic properties 
as shown below. 


3.1. Sources on the symmetry axes 


We consider a source on the axis a = 0. In this case, we 
can find analytic solutions for the lens equation as follows. 
To consider the case of 6 = 0, it is enough to make a 
replacement as a <-> 6 , x <-^ ?/ and 1 -I- e <-> 1 — e. 

For a = 0, Eq. (|^) becomes 



1 - e 


•\/x2 -I- y2 _|_ £.2 


= 0 , 


(31) 


p = 1 -I- ecos2^. 


(27) which means 


In the similar manner to the outer caustic given by Eq. 
(p^), we find 


p = 0 , 


which is the origin in {x,y) (Fig. 2). 
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In the case of a; = 0, Eq. becomes the fourth-order 
polynomial for y as 

E{y) = / - 2by^ + [b^ -h - (1 -f €f]y'^ - 2bc^y + 

= 0. (34) 

Explicit solutions for a fourth-order equation take a 
lengthy form (e.g. van der Waerden 1966). Now, Eq. 
implies 

b-y 


< 0 . 

y 

Using 

E(0) = b'^c^ > 0, 

E(&) = -&2(l + e)2<0, 


(35) 


(36) 

(37) 


we find that E{y) has at least one zero point between 0 
and b. In other words, E{y) = 0 has at most three roots 
for {b-y)/y < 0 . 

Next, we consider the case of Eq. (^. For e = 0, Eq. 
( p 0 | ) means 6 = 0 , so that a ring image appears at 

x^+y‘^ + c^ = l. (38) 

We assume e 7 ^ 0 in the following. Eliminating 

from Eqs. (^0|) and (Is^), we obtain 


y = 


(l-e)b 

2 e ■ 


Substituting this into y in Eq. ([3^, we obtain 
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1 -e 
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3.2. Off-axis sources 

Here, we consider off-axis sources (a 7 ^ 0 and 6 7 ^ 0). In 
this case, Eqs. (^) and ( ^ imply x 7 ^ 0 and y 7 ^: 0. 
Eliminating x"^ + y"^ + from Eqs. (|^ and (|30|), we 
obtain 

(1 — e)bx 


y = 


(43) 


(I -I- e)a — 2ex 

Equation ( p0| ) shows that y is finite for finite 6 , so that 
X 7 ^ (1 + e)al2r from Eq. (|^). Substituting Eq. ( |4^ ) into 
Eq. (|29|), we obtain the sixth-order polynomial for x as 

E{x) = {a — x)^ 

X (x^ -I- c^)[(l -I- e)a — 2ex]^ + (1 — e)^ 6 ^x^ 

— (1 — e)^x^[(l -I- e)o — 2ex]^ 

= 0. (44) 


This equation has at most six real solutions whose analytic 
expressions can not be given by algebraic manners (e.g. 
van der Waerden 1966). As shown below, however, six real 
roots never mean six images. 

In the same manner as for the singular isothermal el¬ 
lipsoid, we obtain a condition for x as 


a — X 
X 


< 0 . 


(45) 


Let us prove that there exists a root between 0 and a. 
We can find out 


F(0) = a*c^{l + ef > 0, (46) 

F(a) = -a^il - e)‘^ < 0. (47) 

Since E{x) is a continuous function, E{x) =0 has at least 
one root in the interval. For (a — x)/x < 0, consequently, 
F(x) = 0 has at most five solutions. Since the polynomial 
is sixth-order, the discriminant is not sufficient to deter¬ 
mine the exact number of roots. Hence, the determination 
is beyond the scope of our paper. 


3.3. Polar coordinates 


Up to this point, we have used 2-dimensional Cartesian 
coordinates: We must solve Eq. (^J) and choose appropri¬ 
ate roots which satisfy the inequality by Eq. (^. Here, 
we adopt 3-dimensional polar coordinates to simplify the 
inequality, as shown below. By taking c as a fictitious third 
dimension, we define 


r = y^x 2 -I- j /2 _|_ (.2 > c, 

T c 
cosw = -, 
r 

X = T sin 41 cos </), 
y = r sin 41 sin (p, 


where we can assume sin4' > 0. Then, Eqs. (29) and 
are rewritten as 
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6 = r 




sin 41 cos 4 >, 
sin 41 sincj). 


(52) 
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We concentrate on off-axis sources (a 7 ^ 0 and 6 7 ^ 0). 
Eliminating sin 41 from Eqs. (p2h and (m), we obtain 


tan (h = 


6 [r- ( 1 -e)] 
a[r- ( 1 -be)]’ 


(54) 


which determines tan^ uniquely for any given r. 
Equations ( |^ ) and (|^) show that r 7 ^ 1 ± e for non¬ 
vanishing a and 5, since x 7 ^ 0 and y ^ 0 mean cos (j> Q 
and sincj) 7 ^ 0. Hence, Eqs. ( p^ ) and ( |5^ ) can be rewritten 
as 


X = 

y = 


ar 

r- (1 - e)’ 
br 

r - (1 -I- e) ’ 
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Substituting these into = a;^ + + c^, we obtain the 

sixth-order equation for r as 

G(r) ^ + 

-a^r^[r - (1 -l- e)]^ - - (1 -I- e)]^ 

= 0, (57) 


which has at most six real solutions. Let us show that there 
are at most five roots compatible with r > c. Using G(c) = 
—a^c^[c — (1 + e)]^ — b‘^c^[c — (1 — e)]^ < 0 and G(— oo) = 
-|-oo > 0 for a continuous function G(r), we find that 
G(r) = 0 has at least one root for r < c. Consequently, it 
has at most five roots for r > c. 

We should note that r and tan (f) are not enough to de¬ 
termine the location of images. A strategy for determining 
the location is as follows: First, we solve Eq. for r > c. 
Next, we substitute r into Eqs. (M) and (^) to obtain the 
image position as {x,y). 
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4. Conclusion 

We have carefully reexamined the lens equation for a cored 
isothermal ellipsoid both in 2-dimensional Cartesian and 
3-dimensional polar coordinates. We have shown that the 
nonlinearly coupled equations are reduced to a single real 
sixth-order polynomial Eq. (jij) or (|^, which coincides 
with the fourth-order equation for a singular isothermal el¬ 
lipsoid as the core radius approaches zero. For the singular 
case, explicit expressions of image positions for sources on 
the symmetry axis are given by Eqs. © and 

©)■ Furthermore, we have presented analytic criteria for 
determining the number of images, which correspond to 
the caustics given by Eqs. (^ and (p^). Consequently, 
analytic expressions for the critical curves are given by 
Eqs. (p^ and (p®. We have shown for the cored case that 
a condition Eq. (|45|) or (^) gives us physical solutions of 
the sixth-order polynomial, which are at most five images. 

The present formulation based on the one-dimensional 
equation (Q or ( ©) enables us to study a cored isother¬ 
mal ellipsoidal lens with considerable efficiency and accu¬ 
racy, in comparison with previous two-dimensional treat¬ 
ments for which there are no well-established numerical 
methods (Press et al. 1988). Particularly for a source close 
to the caustics, the image position is unstable so that care¬ 
ful computations are needed. The amount of computations 
can be reduced by our approach. As a result, it must be 
powerful in rapid and accurate parameter fittings to ob¬ 
servational data. 
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